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Abstract 

We study non-equilibrium steady state transport in scale invariant quantum junc¬ 
tions with focus on the particle and heat fluctuations captured by the two-point current 
correlation functions. We show that the non-linear behavior of the particle current af¬ 
fects both the particle and heat noise. The existence of domains of enhancement and 
reduction of the noise power with respect to the linear regime are observed. The im¬ 
pact of the statistics is explored. We demonstrate that in the scale invariant case the 
bosonic particle noise exceeds the fermionic one in the common domain of heat bath 
parameters. Multi-lead configurations are also investigated and the effect of probe 
terminals on the noise is discussed. 
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1 Introduction 


The study of particle and heat current fluctuations in quantum systems away from equi¬ 
librium attracts much attention both from the theoretical and experimental point of view. 
Such fluctuations generate noise and therefore spoil the signal propagation and detection. 
It is known n-Ei however that current fluctuations carry also useful information, providing 
the experimental basis of noise spectroscopy. The recent progress mu of measurement 
techniques indicates the great importance of this spectroscopy for gaining a deeper insight 
in the mechanism of non-equilibrium quantum transport at the microscopic level. The latest 
conhrmation of this fact comes from the observation of neutral modes [ 6 ] and spin noise 
in fractional quantum Hall edge states via noise spectroscopy. At the theoretical side, the 
derivation of exact results beyond the linear response approximation is fundamental in this 
context, because the current fluctuations are dominated by non-linear effects. 



Figure 1: Two terminal system with scattering matrix §. 

The present paper represents a continuation of our study [ 8 ] of non-equilibrium transport 
at criticality. The focus here is on the particle and heat current fluctuations of systems with 
the structure shown in Fig. [H We consider two {i = 1,2) semi-inhnite leads L* attached at 
inhnity to the heat reservoirs Ri = {f3i,fii). The interaction between the leads is localized 
in the junction and is described by a 2 x 2 unitary scattering matrix S. The system is away 
from equilibrium if the two leads communicate via non-vanishing transmission (off-diagonal) 
elements of S. It is well known that the physical situation, presented in Fig. [T|, is nicely 
described by a Landauer-Biittiker (LB) non-equilibrium steady state constructed from 
the data (/9j,/ii) and §. We have shown in [S] that for scale invariant (critical) §-matrices 
the one-point expectation values of the particle and energy currents in the state can be 
computed in exact and explicit form, which fully takes into account all non-linear effects and 
reveals the presence of an interesting process of transmutation between heat and chemical 
potential energy. Here we pursue further our analysis, deriving the exact two-point current 
expectation values, which allow to determine the particle and heat current fluctuations in 
the state beyond the linear response approximation. In order to test the framework, 
we hrst reproduce and generalize some results concerning the particle noise. In particular, 
we show that this noise depends not only on the parameter /i_ ~ pi — p 2 , but also on 

~ A^i + /^ 2 - The dependence on for /r_|_ = 0 and (3i = (52 has been analyzed in the 
pioneering works of Martin and Landauer [5] and Biittiker (see also nnca). A peculiar 
property of this case is the linear dependence on of the particle current flowing in the 
junction. In order to investigate the impact of the non-linear effects in /i_, we explore in 
this paper the case 7 ^ 0 in a systematic way. This regime is of practical interest, because 
the chemical potentials can be varied and measured directly in experiments. We show in 
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fact that is the parameter which actually controls the noise power in the LB state 
Compared to the case = 0, we observe a relevant noise enhancement for > 0 and 
reduction for < 0. An universal upper bound on the particle noise away from the scale 
invariant regime is also established. Following the main steps of our particle noise analysis, 
we derive and study in detail the less investigated heat noise as well. 

Both fermions and bosons are considered and the influence of the statistics (in particular, 
of the exclusion principle) is discussed. We show that the fermionic noise is a concave function 
of the transmission probability in the junction, whereas the bosonic one is convex. This 
difference provides a new interesting signature for noise spectroscopy. Another remarkable 
manifestation of the statistics concerns the strength of the noise power. We show that for the 
same heat bath parameters in the sate the bosonic noise always exceeds the fermionic 
one. A generalization of the framework and results to junctions with n > 2 leads is also 
presented. In this case we describe the impact of probe terminals on the quantum noise. 

The paper is organized as follows. In section 2 we briefly review the scale invariant 
Schrodinger junction and the concepts of particle and heat current fluctuations. In section 3 
we compute the particle noise power in the fermionic case and describe its general properties, 
focussing on the non-linear effects for 7 ^ 0. Section 4 concerns the heat noise for fermions. 
The bosonic particle noise is derived in section 5. We also perform there a comparison 
between the bosonic and fermionic cases. Multi-terminal conhgurations are discussed in 
section 6 . Section 7 is devoted to our conclusions. In the appendices A-C we address some 
technical issues. 


2 The scale invariant Schrodinger junction 

Referring for the details to [ 8 ], we start by recalling those main features of the Schrodinger 
junction needed in the derivation of the particle and heat noise below. We consider systems 
in which the particle number is conserved (there is no particle production and annihilation) 
and denote by the scattering matrix describing the interaction in the n-particle sec¬ 
tor. It has been argued in [1] that the idea of freely propagating particles along the leads 
accounts remarkably well for the experimental results. In fact, the experiments performed 
with quantum point contacts have shown [T3] that the simple model of independently moving 
electrons, which allow only for the Pauli principle, hts the data with a very good precision. 
In this spirit we assume that = I for n > 1 and keep non-trivial only the one-body 
scattering matrix = §, which drives the system away from equilibrium and is essential 
in this respect. We concentrate mainly on scale invariant (critical) §-matrices because they 
are simple enough to be analyzed explicitly nm and are expected [l5] to incorporate the 
basic universal features of one-dimensional quantum transport phenomena. 

Let us shortly describe now a concrete realization of this general picture. The dynamics 
along the leads, represented by half lines with coordinates {(a;,i),: a; < 0, z = 1, 2} in Fig. 
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( 2 . 1 ) 


[H is fixed by the Schrodinger equatioiJi] 

(i^t + =0. 

The one-body scattering matrix § is determined by requiring that the bulk Hamiltonian 
admits a self-adjoint extension in x = 0. All such extensions are defined [I6]-[TH] by the 
boundary conditions 


+ i(I + ^)ijdx] X, j) = 0 , 


x —>0 


( 2 . 2 ) 


i=i 


where U is a 2 x 2 unitary matrix and A G M is a parameter with dimension of mass. Eq. 
fl2.2p guaranties unitary time evolution and defines a specihc point-like interaction at the 
junction. The explicit form of the corresponding scattering matrix is [T6]-|18j 


§(fc) = 


[A(I - U) - k{l + U)] 

' [A(I - U) + k{I + U)] ’ 


(2.3) 


k being the particle momentum. In order to determine the scale invariant limit of fl2.3p we 
introduce the unitary matrix U diagonalizing U, namely 


UVW = \]A = diag , 


TT TT 

-— < CKj < — . 
2 2 


(2.4) 


From fl2.3p one can easily deduce that U diagonalizes also S(/c) for any k. One hnds 

k + irji k + iri2' 


Sd(^) = lA*S{k)U = diag 


k — ir]i k — iri2 


where 


rji = Atan(aj). 

Scale invariance implies [19] the following alternative 


Vi = 


0 {ai = 0 ), 

oo (tti = 7r/2) , 


Neumann b.c., 
Dirichlet b.c. 


(2.5) 

( 2 . 6 ) 

(2.7) 


Accordingly, the set of scale invariant scattering matrices in the family fl2.3p is given by 

§ = W§dW*, Wef/(2), §d = diag(l,-1) (2.8) 


and the two isolated points § = ±1. The latter are not interesting because there is no 
transmission between the two leads and the system is therefore in equilibrium. For this 

^We adopt the natural units h = c = = \. 
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reason we concentrate in what follows on fl2.8p . The solution of the problem (12.1112.2p is 
given by [ 20 ] 


JL. r°° fjh 

i,(t, = W ■ P.9) 

j=i -fo ^ 

In the fermionic case the operators {ai{k), a*{k) : fc G M, i = 1, 2} generate an anticommu¬ 
tation relation algebra A+ dehned by [ai{k ), aj{p)]+ = [a*{k ), a*(p)]+ = 0 and [ 21 ] 

[ai{k ), a*{p)]+ = 2Tr[6{k - p)6ij + -Fp)], ( 2 . 10 ) 

where * stands for Hermitian conjugation. The deformation fl2.10p of the canonical anti¬ 
commutation relations implements [22]-[23] the interaction in the junction. The algebra 
A- in the bosonic case is obtained by replacing the anticommutators with commutators. 
The explicit construction of the LB steady state as a linear functional over A± was 
given in [20]. This state dehnes the LB representation of A±, which describes the 
non-equilibrium system in Fig. [Tjin terms of {Pi, pi) and §. For convenience we report in 
appendix A the basic correlation functions of {apk), a*{k)} in 

In order to compute the particle and heat noise we hrst recall the form of the particle 
and energy currents 


j^{t,x,i) = ^ - {dx'Pli’] , (2.11) 

j^{t,x,i) = - {dido'll)*) pj - pj* {dtd^'il})]{t,x,i). ( 2 . 12 ) 

The heat current is dehned by the combination 

j^{t,x,i) = j^{t,x,i) - pij^{t,x,i). (2.13) 

The quantum huctuations of any of these currents in the LB state is given by 

Aj^{t,x,i) = j^{t,x,i) - {j^{t,x,i))p^^ , Z = N,E,Q, (2.14) 

where (• • • denotes the fl^^^-expectation value. The noise in the lead Lj can be extracted 
from the two-point function {Aj^{ti, Xi, i)Aj^{t 2 , X 2 , i)) which, because of fl2.14p . equals 
the connected current-current correlator 


{Aj^{ti,xi,i)Aj^{t2,X2,i))fi,i, = (j^(ti,xi,i)j^(t 2 ,a; 2 ,i))““. (2.15) 

Since the energy is conserved in the state the right hand side of fl2.15p depends only 
on the difference fi 2 = — ^ 2 - Accordingly, the noise power at frequency z/ is dehned by the 

Fourier transform 


Pf{P,p,E>;xi,X2;iy) 


dt {j^{t, xi, i)j^ ( 0 , X 2 , i)Yp 


conn 
M ' 


(2.16) 
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Following i-ra. we focus in what follows on the zero frequency noise power 

Pf{l3,ii,S)= (2.17) 

^->• 0 + 

which has the remarkable feature to be a:i^ 2 -independent j2Uj . 

The particle current fl2.1ip satishes at the junction x = 0 the Kirchhoff rule 

j'^(t,0 ,l)+j'^(t,0 ,2 )=0 , (2.18) 

which implies 

F"(/3, S) = Pi'(13, fi, S) = P"(/3, S). (2.19) 

Differently from j^, in general the heat current fl2.13p does not satisfy the Kirchhoff rulej^ 
Therefore, P^ and P® are generally different. Nevertheless, they are simply related. In fact, 
one can obtain P^ from P^ and vice versa by the exchange operation Pi y-)- P 2 and fii fi 2 - 
For this reason we focus below on P®, setting for simplicity 

P«(A/r,§)^Pi'^(/3,/r,S). (2.20) 

In what follows we will hrst derive and study in detail the noises Pf^ and P^ for fermions. 
Afterwards we will summarize the results for bosons and discuss the effect of the statistics. 


3 Fermionic particle noise 


Using the two-point correlation function flA.70p in the representation (see appendix 
A), after some algebra one gets 


1*. t) = + t(1 


(3.21) 


T = |Si 2 p = IS 21 P being the transmission probability and 


P°° rtk k 

(/^, h) = / 73 -+ ^2(fc) - dl{k) - dl{k)\ , 

Jo 

r°° rtk k 

B^{P,ja)= / -- [di{k) + d2ik) - 2di{k)d2{k)] , 

Jo 271 m 

g-/3i [uj{k)-iii] 

diik) =- „, ' —r, z = 1, 2 


(3.22) 

(3.23) 

(3.24) 


is the Fermi distribution in the heat reservoir Pj. The behavior of the particle noise as a 
function of the heat bath parameters {Pi, fip and transmission probability r is fully described 
by fl3.2UI3.24p . We stress that P/^ is well dehned for /9i 2 > 0 and on the whole /ri^ 2 -plane. 


^This is the origin of the energy transmutation in the junction described in [5]. 
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3.1 General properties 


Eq. fl3.2ip provides an useful representation of Indeed, from the explicit form of the 
integrands of Af and one can easily deduce that 

Pf(A/r)>Af(/?,/r)>0, (3.25) 

which, combined with 0 < r < 1, implies that the particle noise is non-negative, 

Pf'^(/5,h,r)>0. (3.26) 

The lower bound (I3.26p is actually a direct consequence of the positivity of the scalar product 
in the Hilbert space "Hlb- 

Concerning the r-dependence, we hrst observe that fl3.25p implies that P/^ is a concave 
function of r. On physical grounds one might be tempted to believe that P/^ is maximal 
at maximal transmission = 1. This is indeed the case only in the range of heat bath 
parameters for which > Bf > Af. If instead B^ > 2Af, the noise P/^ reaches its 
maximum for some < 1. This behavior is illustrated in Fig. HJ 



Figure 2: Plot of the ratio fi, r)/Pf^(/5, fi, 1) for 2Af > B^ > Af (dashed line) 

and Bf > 2Af (continuos line). 


Let us investigate now the dependence of P/^ on the chemical potentials. In order to 
make contact with the previous work i-ia , we set 

/3i = /32 = /3 G 5 ~ 2 ^ ^ ^ ’ (3.27) 

being the average charge density of the two heat baths. The fc-integration in (13.22113.23p 
can be performed exactly (see appendix B) and one hnds 




r f r [l + 2e^(^-+^+) + e2^^-] 

271 f3 1 [1 




P(m_+p, 

(1 — r)(/5/r_) coth(/9/r_) -f (1 — r) coth(/5/i_) In 


+ 


,/3(p_ -M. 

I g/3(p_+P+) 


(3.28) 
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For /i^ = 0 this expression significantly simplifies to 


^f^(/^,h_,0,r) = ^ [r + (1-r)(/?/i_)coth(/?p_)] , (3.29) 

which is precisely the resnlt (with our normalization fl2.1ip of the current) reported in [9]-jl2j. 
For /i_ = 0 one gets instead 


Ff^(/?, 0, /i,, r) =-^. (3.30) 

Because of (13.271) in this case /xi = ^2 and fl3.30p describes therefore the equilibrium current 
fluctuations. In the limit —)■ 0 these fluctuations lead to the Johnson-Nyquist law 

Pf(/3,0,0,T) = ^. (3.31) 

The general expression fl3.28p has a number of interesting properties. First of all is 
an even function of /x_. To our knowledge, the behavior P/^ as a function of /i^ has been 
poorly analyzed previously. It is instructive to fill this gap, because it turns out that the 
noise power depends essentially on the parameter /i^. In fact, one can directly verify that 

> 0, lim P/^(/?,p_,/i ,r) = 0, (3.32) 

hold on the whole parameter space {(3 G M+, /i_ G M, r G [0,1]). Therefore, P/^ is a positive 
monotonically increasing function of the parameter which provides a simple mechanism 
for noise control. In particular, for < 0 (/r^ > 0) the noise P/^ is suppressed (enhanced) 
with respect to the /i^ = 0 value given by (13.291) and studied in [9]-[12]. Fig. [3] displays 
the particle noise for three different values of . The area with green filling represents the 
domain < 0 of noise reduction. A significant suppression of the noise is observed in the 
range /i+ < /i_ < —which is in agreement with the fact that P/^ vanishes in the limit 
—00. 



Figure 3: The noise P/^(l,/r_,/i_^, 1/8) (in units of 1/(3) as a function of fi_ for = 10 
(black dotted line), /i+ = 0 (red line) and /i^ = —10 (blue dashed line). 
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Figure 4: The left panel represents /r_,/i^, 0.9) (in units of 1//3) as a function of /i_ 

for = 20 and (3 = 0.3 (black dotted line), (3 = 0.2 (black dashed line) and (3 = 0.1 (black 
continuous line). The right panel displays P/^(/5,/r_,/i^, 0.9) for /i^ = —20 and (3 = 0.3 
(blue dotted line), (3 = 0.2 (bine dashed line) and /3 = 0.1 (bine continnous line). 


The dependence of P/^ on the temperature 1/(3 is illustrated in Fig. 01 As expected the 
noise increases with increasing the temperatnre both for > 0 (left panel) and /i_^ < 0 
(right panel). 

The shot noise is the zero temperature limit of fl3.28p . Performing this limit, one can 
assnme withont loss of generality that //_ > 0 (eqnivalently pi > /i 2 )- One hnds nnder this 
condition 


ir(l-r)/i_, 

^r(l-r) +/i_), , (3.33) 

0 , /i^ < -/r_ . 

which shows the absence of shot noise [ 22 ] at vanishing transmission (r = 0 ) and at fnll 
transmission (r = 1). We also see that the shot noise 03.331) is a continues function of 
with a discontinnous derivative in = ±/i_. 

Snmmarizing, the general result 03.28P clearly shows that controls the noise power in 
the LB state This feature is of practical relevance because the parameter can be 

directly accessed in experiments. In the next snbsection we will relate the behavior of P/^ 
to the non-linearity of the characteristics of the jnnction for /x^ 7 ^ 0 . 


lim Pf^(/3,/i_,/x+,r) = 

/3^+cxd 


3.2 Particle noise as a function of the current 


In the applications [ 6 ] to noise spectroscopy it is usefnl to express the noise fl3.28p as a 
function of the particle current 


J (/3,/x_,/i+,r) = (j (f,x, 1));3,M = 




1 + e^^^+ 




X _|_ g/3(X‘+-M_) 


r 


(3.34) 


















As already mentioned in the introduction, the value = 0 is very special because in this 
case the current 

Jq = h-, 0, r) = (3.35) 

depends linearly on /i_. Combining this expression with fl3.29p one gets 

271 Jo /t, 0, r) = ^ [t^ + (1 - r)27rl3Jo coth {27r^ Jq/t)] , (3.36) 

which has been tested experimentally in |T3]. The zero temperature limit of f|3.36p gives in 
particular 

hm Pf{(3,271 J^/t,0,t) = {1-t)\J^\. (3.37) 

0 ^ + 00 



0 - 

Figure 5: The current J^(l,/i_, 0.8) (in units of 1/(3) as a function of for /i^ = 2 

(black dotted line), = 0 (red line) and /i^ = —2 (blue dashed line). 


Our goal here is to generalize fl3.36p to the case /i^ 7 ^ 0. It is evident from fl3.34p that for 
7 ^ 0 the current depends non-linearly on , which is illustrated by the blue (dashed) 
and black (dotted) curves in Fig. O It is worth mentioning that such non-linear dependence 
has been experimentally observed in the double-barrier tunnel junctions studied in [28]. In 
fact, the curve [28] relative to these devises resembles very much the blue dashed 

line in Fig. (5] Analogous non-linear behavior has been detected in superconducting tunnel 
junctions in [ 3 ]. 

Solving fl3.34p for /i_ one gets 


/i.(J-) = -ln 


1 _|_ y/l {q2tt0J^/t _ g2/3/i_^+27r/3 J^/t 




(3.38) 


By means of fl3.38p one can eliminate in favor of in the general formula fl3.28p and 
obtain the noise power P/^{(3, fi_{J^), The dependence of this function on is 

displayed in Fig. [Hi The red curve corresponds to /i^ = 0 and describes the behavior 
of fl3.36p studied experimentally for different temperatures and conductances in [13]. The 
new black (left panel) and blue (right panel) curves are obtained for > 0 and < 0 
respectively. We see that for the same value of the current flowing in the junction, 
negative (positive) values of reduce (enhance) the noise with respect to = 0 . 
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Figure 6: P/^ (in units of 1//3) as a function of for (3 = 1 and r = 1/6, represented for 
= 0 by a red line in both panels. The same function for /i^ = 1 (black dashed line) and 
= 3 (black dotted line) in the left panel, and = — 1 (blue bashed line) and = —10 
(blue dotted line) in the right panel. 

Let us observe that according to (I3.38p the value = 0 implies /i_ = 0 and leads 
therefore to the equilibrium current fluctuations given by fl3.30p . 


3.3 Relaxing the condition of scale invariance 

We conclude our study of the particle noise by considering momentum dependent scatter¬ 
ing matrices. Scale invariance in the junction is broken in this case and in general the 
fc-integration cannot be performed in exact and explicit form. Nevertheless, one can es¬ 
tablish some useful estimates. Let us introduce for this purpose the maximal and minimal 
transmission probabilities 

Tmin = min{|§i2(fc)P : k> 0 }, r^ax = max{|§i2(fc)P : k> 0 }, ( 3 . 39 ) 

which exist because 0 < |Si 2 (fc)P < 1 by unitarity. Then P/^ satishes the following inequal¬ 
ities 


Pf^(/3,/i,§) > r2.^klf(/3,/i) + r„ia(l-w)Pf(/9,/i), (3.40) 

Pf{(3,fI,T) < r^,,,,klf(/3,/i) + rmax(l-^min)5f(/3,/i), (3.41) 

where Af and Bf are the ^-independent integrals (13.22113.231) . The estimates (13.40113.411) 
follow only from unitarity and in this sense are universal. 

Let us consider at this point the scattering matrices (12.3p generated by the boundary 
condition (12.21) . Inserting in (12.5p the most general 2x2 unitary matrix U, one obtains [20] 
the transmission probability 

^ (3.42) 
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where r]i are given by fl2.6p and the angle 9 is the only parameter which remains from the 
four ones characterizing a generic U G U{2). From fl3.42l) one gets 


^min 


0 , 




iVi -r 72 )^sin^( 6 >) 

(Ihil + \V2\y 


(3.43) 


The combination of (I3.4n|l3.4ip with (I3.43p implies 

0 < /i, §) < 9) p) + 9) Sf (/3, p), (3.44) 

which conhrms (I3.26p away from criticality and provides an upper bound on the particle 
noise for the whole family (12.3p of fc-dependent scattering matrices. 


4 Fermionic heat noise 

The above analysis can be extended to the heat current fl2.13p . One finds 

F/^(/3,/i,r) = r2ylf(/3,/i) + r(l -r)Ff^(/3,/i), (4.45) 

where 

^ dk k 

Af{l3,fi)= - - [u{k)-iJ,i\ [di{k) + d2ik) - dl{k) - dl{k)] , (4.46) 

Jo 

^ dk k 

B^{/3,iJ,)= / -- [u{k) - Hi] [di{k) + d2{k) - 2di{k)d2{k)] . (4.47) 

Jo 

Analogously to the particle noise, the integrals fl4.46114.47p satisfy 

(4.48) 

implying that the heat noise is non-negative, 

F«(/3,/x,r)>0. (4.49) 


The r-dependence of is similar to that of the particle noise. In fact, F® is maximal 
at = 1 if 2Af > B^ > Af. If instead F® > 2Af, the heat noise P^ reaches its maximum 
for < 1. The dependence of the heat noise fl4.45l) on fd and /i^ is much more involved. 
For the sake of conciseness we report the explicit form of P^ in appendix C, summarizing 
here its basic properties. One can deduce from flC.87D that 

lim Ff^(/3,/i_,/i^,r) = 0, (4.50) 

which implies that F® is an increasing function of h+- Differently from P^ however, F® is 
not symmetric in jj,-- 
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Figure 7: The noise /i_,/i^, 1/2) (in nnits of 1//3) as a fnnction of //_ for /i^ = 20 

(black dotted line), /i^ = 0 (red line) and /i^ = —20 (blue dashed line). 


The plots in Fig. [7] give an idea about the behavior of for different values of /i^. The 
heat noise power is reduced for negative and enhanced for positive exactly like the 
particle noise P/^. 

The zero temperature limit of flC.87p gives (for /i_ >0): 


lim P®(/3,/x_,/i^,r) = <^ ^r(l-r) (/x+ +/i_)^ 

5—>-00 I ^ 

lo> 


> h- , 

-/i_ < < /i_ , (4.51) 

< —/i_ . 


At equilibrium (/x_ = 0) one hnds 

r{;3Xe'^'‘.-2(l + e-^>-qp^nhl+e'^'‘b+L.q-e'’'-p]} 

t', ' ^ ^ 

Setting in addition /i^ = 0 one has 

P«(/?,0,0,r) = ^, (4.53) 

which is the counterpart of the Johnson-Nyquist law concerning the heat noise. 


5 Bosonic particle fluctuations 


In the bosonic case one gets from fl2.15p - fl2.17p the particle noise 

P/"(/3, /X, r) = h) + r(l - r)Pb^(/3, /x), 

where now 

r°° rtk k 

/ X- [di{k) + d2{k) + dl{k) + dl{k)~\ , 

./n ^71" rn, 


(5.54) 

(5.55) 
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r°° dk k 

/ 7. - [di{k) + d2{k) + 2di{k)d2{ky\ , 

and di{k) is the Bose distribution 


(5.56) 


di{k) 


1 — Q-9i[^^{k)-^ii\ ’ 


i = 1,2. 


(5.57) 


A new characteristic feature of the bosonic case is the presence of singularities in the in¬ 
tegrands of fl5.55115.56p aX k = In order to avoid them and deal with meaningful 

expressions, we impose 

/xi < 0 and /i 2 < 0 /i^ < 0 and /i+ < /i_ < —/i+ . (5.58) 


This assumption shifts the singularities away from the range of integration. We stress that all 
our results about bosonic systems hold only in the open cone C = < 0, /i_^ < /i_ < —d+} 

defined by fl5.58p . The integrals fl5.55115.56p are well dehned there and satisfy 

A^(/3,/i)>i?b^(/?,/i)>0. (5.59) 

The bosonic noise power is therefore non-negative like the fermionic one. However, 
differently from the fermionic case, P-^ is a convex function (see Fig. [8]) of the transmission 
probability r. This fact follows directly from fl5.54115.59p and provides a possible tool for 
detecting the statistics in noise spectroscopy. 



Figure 8: Plot of the ratio = P^ {(3 , fi, t) /P^ {(3 , fi, 1). 


The cone C represents the common domain where both Pf^ and P^ are well defined and 
can be therefore compared. From the inequalities flB.77IIB.78p in appendix B it follows that 
for the same heat bath parameters the bosonic noise always exceeds the fermionic one, 

P^{(3,f,)>P,^{(3,ix), (/i+,/r_)eC. (5.60) 

This remarkable feature illustrates the deep role [T] of the exclusion principle in reducing 
the quantum noise in the domain C, where the bosonic transport is well dehned. The result 
(15.601) can hnd interesting applications in quantum transport. 
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For (3i = (32 = 13 the explicit expression of in C is (see appendix B) 




r 

27r/3 


r [2e^(M_+/^+) _e2/3M_ _ 

I^X — e^(/^-+A‘+)j l^x _ g/3(M_-/i+)j 


(1 — t){(3^_) coth(/3p_) — (1 — r) coth(/5/i_) In 


X — g/3(/^_+/^+) 

1 

g/^M- _ g/3A‘+ 

i 


(5.61) 


The presence of singularities outside of the open cone C is manifest in (I5.6ip . which diverges 
on the boundary /r_ = and on the tip = 0 of C. The lines /i_ = 0 with 7 ^ 0 
belong to C and 

Ph (/^5 0,/r+,r) = _e/3M+) ’ 

which describes the equilibrium bosonic particle fluctuations. 



y- 

Figure 9: The noise Pi(^(l,/i_,/r^, 1/2) (in units of 3/(3) as a function of /r_ for = —1 
(red line), = —1.5 (dashed blue line) and /i^ = —2 (dotted blue line). 


Like for fermions, is a symmetric function of /i_, which increases monotonically with 
G (— 00 , 0) and vanishes in the limit — 00 . Fig. |9] illustrates the behavior of P^ 

in the variable /i_ for three different values of /i,^. With the substitution Af 1 —)■ and 
I—)■ the bounds (13.40113.41113.44p hold for non-critical bosonic junctions as well. 


6 Multi-terminal junctions 


We generalize here the above results to junctions with n > 2 leads, concentrating on the 
particle noise power for fermions. Plugging flA.TOp in fl2.16p one hnds in the i-th terminal 


P(''(/?,A,S;i) = 

__ J (1 _ 2\^u\^)d,ik) [1 - d,ik)] + \^ii\"diik) - 

71 m \ 




1=1 


(6.63) 
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In order to illustrate the physics behind this complicated formula, which depends on many 
parameters, it is instructive to consider the case n = 3. One can interpret the lead i = 3 as 
a probe terminal and study for instance its influence on the noise in the lead i = 1. Setting 
T = |Si 2 p and a = |§i 3 p one gets from fl6.63p 


/i, §; 1) = (/3, /r; r)(i, 2 ) + P/'' (/?, /x; a)(i,3) + P/^ (/?, hi P a)(i,2,3), 


(6.64) 


where P/^(/5,/i; r)(i_ 2 ) coincides with the two-lead expression fl3.2ip . P/^(/3,/i; (t)(i^ 3 ) is still 
(13.211) but with (32 '—> /^s, Ai 2 '—> '—> o' and Pf^(/5, p; r, cr)(i, 2 , 3 ) is a term mixing all 

three reservoirs with the form 


r°° dk k 

Pf^(/?,/r;r,cr)(i, 2 , 3 ) =-2rcr / -- [di{k) - d 2 {k)][di{k) - d^ik)] . (6.65) 

Jo m 

As expected, fl6.64p is symmetric under the exchange 2 t-)- 3. For (3i = (32 = (J^ = (3 the 
fc-integrals can be computed explicitly. One has 


Pf^(/3,/i,§;l) = Pf^(/3,p_,/i^,r) + 

^f^(/3,h+ +h- -h3,h+ +h- +h3;o-) + Pf^(/?,//+ +h_,h+ -h-,h3;pcr), (6.66) 


where the hrst two terms are given by fl3.28p and 

Pf^(/3,Pi,P2,h3;pc^) = 


rcr 

TT/d 


In (l -|- e^^^) e^'^^ In (l -|- — e^^^ In (l -|- 


X -|- Q| 3 ^J.l g/3Mi — QPfJ-2 

,/3/ii In ^x -|- e^^^) — e^^® In (l -|- e^^^) e^^^ In (l -|- e^^®) — e^^® In (l 4- e^^^) 


3(3^1 — p/3m3 


=,/3At2 _ p/3M3 


(6.67) 



Figure 10: The noise P/^(l,/i_,= 0,/r3 = 10, r = 0.4, a) (in units of 1/(3) as a function 
of /r_ for cr = 0 (red line) and a = 0.5 (dashed blue line). 

If one isolates the probe terminal by setting a = 0, eq. fl6.66p exactly reproduces the 
noise fl3.28p of the two-terminal junction. This feature is illustrated by the continuous (red) 
curve in Fig. [TOl For a ^ 0 instead, the probe terminal affects the noise as shown by the 
dashed (blue) line in the same hgure. We see that turning on a probe terminal reduces the 
noise outside of some neighborhood of fi_ = 0. The addition of probe terminals represents 
therefore a viable mechanism for decreasing the noise. 
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7 Outlook and conclusions 


We studied in this paper some non-linear effects in the particle and heat current fluctuations 
in scale invariant quantum junctions in the Landauer-Biittiker non-equilibrium steady state. 
These effects are associated to the non-linear dependence of the particle current on 
the chemical potential difference fi_. In the case with two heat baths our findings can be 
summarized as follows: 

(i) Depending on the value of the sum /i^ of the chemical potentials, the non-linear 
effects lead to reduction < 0) or enhancement (/r^ > 0) of both particle and heat noise 
for fermions; 

(ii) In the bosonic case the quantum transport is well defined for /i^ < 0 and decreases 
with decreasing of /i^. In the common domain < 0 and with the same values for the heat 
bath parameters the bosonic particle noise always exceeds the fermionic one; 

(iii) Another difference between fermions and bosons concerns the behavior of the particle 
noise as a function of the transmission probability in the junction. For fermions this function 
is concave, being instead convex for bosons. 

In this context we established an universal upper bound on the particle fluctuations away 
of criticality. We considered also multi-terminal junctions, investigating how the chemical 
potential of a probe terminal affects the particle noise. 

The above results have been extracted from the two-point current correlation functions. 
The theoretical progress in the study of mesoscopic systems has led to the concept [2U]-PT] 
of full counting statistics, in which a more detailed information about the non-equilibrium 
particle transfer is obtained from the n-point current correlators with n > 2. We expect that 
the extension of our work to this case will shed additional light on the non-linear character 
of the transport in quantum junctions. 


A The Landauer-Biittiker representations 


For a detailed construction of the LB representation of the algebra A± we refer to [20] . 
The two-point function {a*{p)ai{k)) in is given by 


{a*{p)ai{k))= 27rd{k - p) 


e{k)6,,df{k) + 9{-k)J2^u d^{-k) 


1=1 


+27rS{k + p) [6{k)df{k)S*j + 6{—k)Sijdf{—k)] , 


(A.68) 


were d'l{k) is the Fermi distribution fl3.24p and d^ {k) is the Bose distribution 05.571) . The 
explicit form of {ai{k)a*{p)) is obtained from 0A.68P by the substitution 

df{.k) I —^ 1 T df{k) . (A.69) 
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As well known [32] , employing the relation fl2.10p one can express a generic n-point correlation 


function as a polynomial of the two-point correlators 
Using this property and flA.hSp one finds 


{a*Ap)ai{k ))and {ai{k)a*Ap)) 


M- 


Xi, X2, * 2 ))““ = 

^ 4{h)[lTdt{kP] 

l,m=l 

x{xH,{h]Xi) [d^Xhm] {k2;xi) - [d^Xu,] ih] Xi)xhm{h] Xi)^ 

x{x*mi2if^2]X2) [d:cXi2l] {h] X 2 ) - [dxX*mi^] [h] X2)Xi2l{kl] X2)^ , (A.70) 

where y is the matrix 

y(fc;x) = + (A.71) 

The correlation function xi, a^ 2 A 2 ))^°“ is obtained from flA.70p by inserting 

the factor [uj{ki) + 02 (^ 2 )]^ in the integrand. 



B Particle noise integrals 


Adopting the variables 


y = 


r = A 1 /A 2 , 


Oj = e 


i = 1,2, 


the integrals 



and (15.55115.56P can be written in the form: 





aij/^ ^ 0.2 

_(?/’' +ai )2 ^ ( 2 /+ 02)2 
ai + a2y''~^ 
iy'' + ai){y + a 2 ) ’ 

' , Q 2 
_(?/’'-ai )2 ( 2 /- 02)2 

ai + a2y''~^ 

( 2 /'’ - ai){y - 02 ) ■ 






(B.72) 

(B.73) 

(B.74) 

(B.75) 

(B.76) 


These integrals are well dehned for r > 0 and /ij < 0 (oj > 1). Moreover, in this domain 


/lb /If 


1 /•' r 2 a?|/ 2 '-i 2aly 

J_ f\ 2/(Qi + a2y''~^)‘^ ^ p 

^^ 2 ], \y^^-aiW-al)- ’ 


> 0 , 


since the integrands in 



are non-negative. 


(B.77) 

(B.78) 
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Performing the integration in flB.73IIB.75j) one finds 




1 1 
+ 


27r/?2 + 1) (02 + 1) 


+ 


2vr/52 [r(ai - 1) (02 - 1)J 


(B.79) 

(B.80) 


The explicit form of flB.74IIB.76p for generic r is not known. For r = 1 (/3i = (32 = fi) one 
has 


Bf 


27 rp (ai — 02) 

1 (oi + 02 ) 

—n 7-7 

27rp (ai — 02 ) 


CI 1 CI 2 “1“ CL\ 
(Xi(l2 + CL 2 
Oj\Oj2 — ^2 
Qj\(l2 — Ci\ 


(B. 8 I) 

(B.82) 


C Heat noise integrals 


The integrals fl4.46114.47p can be computed explicitly for = (32 = (3. Using the variables 
flB.721) one hnds 


= 

^ 277/^3 


1 


[C(ai;pi) + C'(a2;pi)] , 


where 


C'(a;pi) = +2/jpiln 

a + 1 


a 

— 2Li2 

1' 

Cl “h 1 

a 



T>(a;pi) = (/?pi)^ln 


2/3piLi2 


+ 2Li3 


Cl T 1 

Substituting these expressions in fl4.45p and restoring the variables p± one gets 

^2 


(C.83) 

(C.84) 

(C.85) 

(C.86) 




’'(1 - t) 

277/93 


g-9(M++M_) _j_ g-/3(7t+-M_) 
g-4(M++M_) _ g-/3(M+-M_) 


[D (e-^('^+ ); p^ + p_) + D (e-^W ^; p+ + p_) ] . 

(C.87) 
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